3D 打 印 进程 中 曲率 对 直立 墙 结构 力学 的 影响 


刘 轩 廷 19 孙 博 华 *D 


!( 西 安 建筑 科技 大 学 理学 院 , 西安 710055) 
西安 建筑 科技 大学 土木 工程 学 院 , 西安 710055) 
i 西安 建筑 科技 大 学 力学 技术 研究 院 , 西安 710055) 


摘要 对 于 打印 墙 结构 往往 因为 稳定 性 而 失效 的 问题 ,本 文 提出 了 采用 直立 波纹 墙 的 方法 , 即 适当 增加 墙 体 治 轮 
廓 方向 的 曲率 ,以 求 提 高 打印 墙 稳定 性 .本 文通 过 将 波纹 墙 的 单个 波纹 视 为 单 胞 进行 研究 , 取 曲 率 为 常数 ,并 进行 
结构 闭合 处 理 ,以 求 稳定 性 最 大 化 , 即 分析 柱 壳 形 直立 墙 在 3D 打 印 进 程 中 的 力学 性 能 .此 模型 利用 Goldenveizer- 
Novozhilov 壳 体 理 论 ,加 入 打印 进程 参数 ,包括 打印 速率 、 打 印 材料 固化 特征 、 柱 壳 几 何 特征 、 以 及 自重 作 

影响 ,对 圆柱 壳 的 两 种 破坏 机 理 : 弹 性 屈曲 和 塑性 破坏 进行 分 析 . 利 用 Suiker22l 首次 提出 的 无 量 纲 参数 入 ( 暂 
称 Suiker 数 ), 对 柱 壳 形 直立 墙 弹性 屈曲 与 塑性 破坏 间 的 竞争 关系 进行 描述 .将 结果 与 文献 5 进行 比较 ,发 现 当 参 
数 完全 相同 情况 下 ,将 柱 壳 形 直立 墙 比 矩 形 直 墙 结果 对 比 发 现 ,稳定 性 提升 了 两 们 有余. 同时 给 出 了 通过 曲率 的 改 
变 达到 此 类 参数 下 柱 壳 结 构 稳 定性 最 大 化 的 方法 , 即 通过 拟 合 的 手段 ,得 到 失效 高 度 随 曲率 变化 函数 , 据 此 对 失效 
机 制 区 域 进行 划分 .这 表明 该 模型 可 用 于 探讨 打印 进程 中 曲率 对 直 墙 结构 的 力学 性 能 影响 ,并 可 通过 调整 参数 对 
结构 作出 优化 .同时 ,通过 引入 曲率 ,将 直立 墙 治 轮廓 方向 打印 成 波纹 墙 , 对 于 提高 打印 墙 结构 稳定 性 是 一 个 简单 
可 行 的 方案 . 
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Analysis of mechanical performances of cylindrical wall in 3D printing processes 
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Abstract In order to enhance the overall structural stability of 3D printed concrete wall, this paper propose a novel 
scheme to produce a wavy wall with curvature along its contour direction. To validate the idea of the scheme, a single 
wave wall is set as a cylindrical wall and its mechanical performances are analyzed. The mathematical 3D printing model 
of the wavy wall is formulated by the Shell theory while taking into account of additional parameters of the printing 
processes, the model will be used to analyze the two failure mechanisms of the cylindrical wall: elastic buckling and 
plastic collapse. Compared with the results of Suiker’s straight wall, it is found that when the parameters were same, the 
stability of the cylindrical wall is more than twice of the rectangular wall. Our studies indicate that it’s a feasible scheme 


to improve the printed structural stability by increasing curvature. 
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引 


了 中 


3D 打 印 的 概念 是 由 Charles Hull 于 20 世 纪 80 年 
代 初 提出 ,其 也 可 以 称 为 增 材 制造 AM), 快 速成 
型 RP), 无 国体 成 型 (SFP) 等 .自问 世 以 来 至 今 30 余 年 
间 ,3D 打 印 技术 发 展 尤为 迅速 ,并 取得 了 巨大 的 成 
功 .其 应 用 包括 生物 医学 工程 ,车 辆 工程 ,机 械 工程 , 航 
空 航天 工程 ,食品 工程 ,建筑 土木 工程 等 [77 术语 ” 增 
材 制造 "最 终 由 美国 材料 协会 (ASTMD) 选 择 , 以 与 传 
统 ” 减 材 制造 " 相 区 分 ,并 于 2013 年 将 增 材 制造 (AM) 技 
术 分 为 七 类 :(1) 材 料 挤 出 、(2) 粉末 床 融 合 、(3) 还 原 


打印 层 数 ,并 为 找寻 最 佳 打 印 参数 集 给 予 了 理论 指 
导 . 但 从 其 结果 发 现 , 矩 形 布置 的 墙 体 结构 稳定 性 仍 
有 不 足 , 即 发 生 结构 届 曲 而 失效 ,而 不 是 达到 材料 塑 
性 极限 发 生 破坏 .由 于 稳定 性 研究 对 于 提升 材料 利用 
率 ,快速 成 型 和 增强 结构 的 承载 力 等 方面 至 关 重 要 。 
特别 是 ， 在 不 改变 打印 速度 的 情况 下 ,使 用 混凝土 
等 流动 性 强 固化 慢 的 3D 打印 材料 时 , 如何 提高 打印 
墙 结构 的 整体 稳定 性 是 个 关键 科学 问题 。 

为 了 解决 这 个 问题 ， 我 们 提出 了 一 个 简单 可 行 
的 普 适 方案 ， 即 适当 增加 直立 墙 沿 轮廓 方向 的 曲率 
将 直立 墙 沿 轮廓 方向 打印 成 波纹 形状 ， 打 印 方案 见 


光 聚 合 、(4) 材料 喷射 、(5) 粘 结 剂 喷射 、(6) 片 状 分 层 
和 、(7) 定向 能 量 沉积 鲜 . 其 主要 形式 是 通过 建立 三 维 
模型 数字 文件 ,利用 软件 将 此 模型 进行 切 分 (精度 取 
决 于 打印 分 状 率 ), 最 后 使 用 打印 设备 进行 逐 层 打印 
成 型 中 AM 技 术 相 对 于 传统 制造 的 优点 是 :能 制造 复 
杂 几 何 形状 的 部 件 ,减少 材料 浪费 ,减少 能 量 损耗 , 增 
加 面 对 客户 的 灵活 性 等 00-11. 

Bourell 等 65 人 (4 于 2009 年 开展 了 关于 AM 技 术 
发 展 规 划 研 讨 会 ,探讨 了 关于 AM 技术 发 展 的 重要 
方面 ,包括 :设计 ,过 程 建 模 与 控制 ,材料 、 打 印 进程 
与 打印 机 器 ,生物 医学 应 用 ,能 源 及 可 持续 发 展 .关于 
增 材 制造 (AM) 的 多 元 化 发 展 ,Behrokh Khoshnevis 提 
出 “轮廓 工艺 ”, 即 混凝土 3D 打印 技术 .此 技术 由 
于 材料 的 局 限 性 ,经 济 性 等 原因 ,尚未 挖掘 出 其 真正 
潜力 ,但 是 在 某 些 极端 条 件 下 ,有 效 的 弥补 了 传统 建 
筑 的 弱项 ,如 在 防疫 方 舱 医院 的 应 用 ,于 其 他 星球 建 
造 栖息 地 等 (3-191， 目 前 就 3D 打印 进程 中 的 力学 性 
能 分 析 相 对 较 少 ,国内 外 主要 研究 集中 在 材料 科学 
方面 07-20. 但 由 于 3D 打印 进程 中 所 涉及 参数 的 复杂 
性 与 多 样 性 ,包括 打印 物体 的 几何 特征 ,时 间 效 应 对 
材料 强度 的 影响 ,不 均匀 的 自重 载荷 ,喷头 的 移动 速 
率 ,不 可 避免 的 几何 缺陷 ,通常 来 说 ,打印 参数 的 设 定 
采用 试 错 的 方式 进行 ,而 针对 混凝土 构件 在 打印 进程 
中 的 力学 性 能 分 析 的 研究 更 是 少 之 又 少 . 

根据 A.S.J. Suiker22 针 对 混凝土 ?D 打 印 进程 中 
的 力学 性 能 研究 发 现 ,弹性 屈曲 与 塑性 破坏 是 结构 
破坏 的 主要 机 制 户 -2.Suiker 针 对 直 壁 结构 3D 打 印 进 
程 中 的 力学 性 能 进行 分 析 , 并 首次 给 出 了 相应 的 力学 
模型 ,在 引入 混凝土 固化 性 能 与 打印 参数 下 ,将 参数 
模型 结果 与 试验 数据 进行 对 比分 析 , 发 现 可 较 好 地 
分 析 结 构 在 打印 进程 中 的 破坏 行为 ,预测 直 壁 结构 


图 1 。 直立 波纹 墙 每 个 波纹 , 其 可 以 是 圆柱 形 也 可 以 
不 是 ， 同 时 其 曲率 也 可 以 变化 。 为 了 验证 我 们 的 广 
案 , 本文 对 于 柱 壳 形 直立 墙 在 3D 打 印 过 程 中 力学 性 
能 进行 了 研究 。 


图 1: wavy-crinkle-crankle-wall 演 化 示意 图 . 


具体 讲 , 在 第 2 节 , 推 导 柱 壳 的 压 曲 控制 方程 .在 
第 3 节 , 描 述 了 3D 打印 进程 中 柱 壳 结 构 可 能 发 生 的 两 
种 失效 机 制 .第 4 节 , 利 用 数值 计算 方法 得 到 了 两 种 固 
化 参数 影响 下 的 屈曲 结果 ,并 使 用 Suiker 数 对 柱 壳 墙 
的 两 种 失效 机 制 间 竞 争 关系 进行 评判 .第 5 节 , 采 用 与 
文献 的 完全 相同 的 打印 参数 ,对 由 曲率 控制 的 柱 壳 结 
构 稳定 性 和 由 材料 属性 控制 的 塑性 极限 进行 了 对 比 
分 析 . 


如 


1 非 均 匀 内 力作 用 下 的 非 均 质 柱 壳 


如 图 2 所 示 , 我 们 取 波 纹 墙 的 一 个 单 胞 进行 分 
析 , 其 曲率 是 可 以 变化 的 , 即 上 = Km, r = variable. 但 
当 我 们 取 曲 率 为 常数 , 即 r = constant ,并 进行 结构 闭 
合 处 理 ,明显 将 得 到 一 个 结构 稳定 性 与 承载 力 最 大 化 
的 柱 壳 形 直立 墙 . 
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图 2: 直立 墙 演 化 示意 图 . 


在 3D 打 印 进 程 中 , 柱 壳 形 直立 墙 会 因为 材料 固 
化 行为 影响 ,在 垂直 方向 呈 非 均 质 性 .此 外 , 柱 壳 结构 
承受 的 载荷 为 自重 载荷 ,其 应 力 在 垂直 方向 (自重 方 
向 ) 同 样 是 不 均匀 的 .与 之 对 应 的 是 , 需 基于 平衡 方程 
与 边界 条 件 ,对 此 非 均 匀 内 力作 用 下 的 非 均 质 柱 壳 的 


力学 模型 进行 描述 .因此 ,本 节 通 过 柱 壳 的 势能 最 小 
化 原理 将 其 压 届 微分 方程 进行 推导 . 
1.1 势能 
图 3: 柱 壳 在 自重 作用 下 示意 图 . 
如 图 3 所 示 , 柱 壳 半 径 为 R ,高 度 为 1 ,均匀 厚度 有 h 
. 轴 对 称 情况 下 , u(x,0) = u(x),v = 0, w(x,9) = w(x) . 根 


据 Goldenveizer-Novozhilov 壳 体 理论 ,应 变 位 移 关系 可 
以 表示 为 : 
Er=Ux, 82y= 


(1) 
Xx = Wxx; Xy 二 0, Xxy = 0. 


其 中 = 倚 , wx = 鱼 , x 为 纵向 坐标 , y 为 环 向 坐标 ， 
By Ey, Ey 为 壳 体 中 面 线 | 性 应 变 分 量 , yw 好 ,Xe 为 
率 变 化 ,u,v, w 分 别 为 纵向 , 切 向 , 径 向 位 移 . 

柱 沉 总 应 变 能 UU ,由 薄膜 应 变 能 UU, ,与 中 面 拉 伸 
弯曲 应 变 能 UU 组 成 : 


U = Us;+ U;, (2) 


2 


y 
2 十 Eg + 2yexEy | dydx, 


(3) 


,2 


Up = {J > ps + 好 +2yxxXy + 2(1 一 wxs| dydx, 
(4) 
星 号 ”*” 代 表 材 料 参 数 在 空间 中 是 非 均 匀 的 (为 方便 
起 见 , 取 v 为 常数 )，E, 为 非 均匀 弹性 模 量 , 非 均匀 弯 
曲 刚度 由 下 式 给 出 : 


| hy 


~ 120 
柱 壳 届 曲 过 程 中 自重 所 做 的 功 可 表示 为 : 


W = ;| | psn DPF + (RF lads (6) 


合 公 式 (2)(6), 可 以 得 到 总 势能 I : 


-fs 


+ 上 本 e + 好 +2yxxXy + 2(1— wxs| dydx 
1 2 2 
一 由 SA 一 由 ba 十 wa| dydx. 


1.2 平衡 方程 与 边界 条 件 


将 公式 (7) 进 行 势能 最 小 化 处 理 ,可 以 得 到 非 均匀 
内 力作 用 下 的 非 均 质 柱 壳 的 平衡 方程 : 


1 
oIl = 3 上 De[w abw xx + WyyOW,yy 


2(1 一 Y 
R2 


13 | { te a 
5 Tt 了 Wow + RWOLx 


2(1 一 
A 
RR R? 


1 
号 上 小 pgh(x 一 D[vrovr+wWwxrow rdydx， 


妆 
Er + Ey 十 27ExEy 


y 
7 了 | dydx 


(7) 


) 
+ yw,xxOW,yy + YW,yyOW,xx + WryOw,xy Jdydx 


uyou,.y]Jdydx 


(8) 


当 处 理 柱 壳 轴 对 称 届 曲 问题 时 .将 (D 式 带 入 ,进行 位 
移 变 分 可 得 柱 壳 压 届 微 分 方程 : 


she — Dw.r) x = 0, 


1 
(CD 页 [oa 
(9) 
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此 外 , 柱 过 沿 壳 边缘 x = 0 和 x=1 的 边界 条 件 分 别 为 


wlx=0 = 0， Wixlx=0 = 0， 
(10) 


Wxxlx=1 三 0， (DiWw,xx),xlx=1 =0. 


2 ”3D 打印 进程 中 柱 壳 的 失效 


plastic collapse 


! 


elastic buckling 


4: 失效 机 制 示 意图 . 


如 图 4 所 示 , 在 3D 打 印 进程 中 可 能 发 生 两 种 失 
效 机 制 ,分 别 为 弹性 屈曲 和 塑性 破坏 .其 中 当 材 料 
硬化 速率 比 打印 速率 相对 较 慢 情况 下 ,在 结构 底部 
会 出 现 较 大 变形 ,文献 59 中 试验 观测 ,第 一 层 材料 由 
于 受 约束 变形 较 小 ,而 第 二 层 在 结构 塑性 破坏 情况 
下 ,最 终 厚 度 将 达到 第 一 层 两 倍 以 上 ,可 称 这 种 现 
象 为 "Elephant foot deformation”. 这 种 现象 的 合理 利 
用 ,可 在 一 定 程度 上 增强 结构 稳定 性 

在 本 节 中 建立 了 这 两 种 失效 机 制 的 控制 方程 . 弹 
性 屈曲 的 控制 方程 可 以 由 第 一 节 中 所 导出 的 压 屈 微 
分 方程 和 边界 条 件 所 得 .引入 打印 材料 固化 过 程 和 印 
刷 速度 有 关 的 时 间 效 应 的 影响 因素 .本 节 最 后 给 出 了 
弹性 届 曲 与 塑性 破坏 的 竞争 标准 . 


2.1 打印 时 间 因 素 对 柱 壳 的 屈曲 影响 


图 5: 柱 壳 打印 过 程 示意 图 
,底部 为 拉 格 朗 日 坐标 系 x, 打印 喷头 处 为 欧 拉 坐 标 
系 X. 


由 于 柱 壳 体 长 度 ! 随 打印 进程 增长 ,其 增长 简化 
为 整个 柱 壳 结构 沿 x 轴 方 向 按 恒定 速度 多 = i 的 连 
续 增 长 过 程 .其 中 速度 i 由 打印 进程 中 的 参数 决定 , 它 
可 以 从 流量 关系 8 = mp7 得 到 , 即 


;__0 
i= 11 
AT (11) 


其 中 为 打印 材料 从 喷头 每 单位 时 间 排 出 体积 ，w 
为 喷头 水 平移 动 速 度 , Ti 为 打印 单个 材料 层 所 需 时 
间 , h 为 壳 体 厚度 .如 图 5 所 示 , 打 印 过程 可 以 通过 采用 
附着 在 喷嘴 末端 的 欧 拉 坐标 系 进行 判定 .其 中 欧 拉 坐 
标 X 与 拉 格 朗 日 坐标 x 的 关系 : 1 = x+X, 可 以 得 到 
欧 拉 坐标 X 如 下 : 


X= (x,t) =x-), (12) 


一 般 情况 下 可 以 认为 打印 速率 是 固定 的 , 即 i = 
constant, 并 设 时 间 1(0) = 0, 有 关系 1= it. 这 样 上 式 可 
以 简化 成 

X= R(x,t) =x— i, (13) 


其 中 上 代表 时 间 . 如 果 打印 速率 是 变化 的 问题 就 比较 
复杂 了 ,将 在 今后 的 研究 中 考虑 . 

考虑 材料 在 点 x = 0 处 ,固化 过 程 中 弹性 刚度 随 
时 间 演 化 可 表示 为 : 


E(x = 0,7) = $&,(1) Eo, (14) 


其 中 8 = 8.07) 为 固化 特征 函数 , Eo 为 材料 从 喷嘴 喷 
出 的 初始 刚度 .如 文献 中 所 示 , 可 以 得 到 欧 拉 坐 标 系 


下 目前 最 具 代表 性 的 两 种 固化 特征 函数 。 
(a) 线性 固化 特征 函数 


亲本 =1- 台 (15) 


(b) 指数 衰减 固化 特征 函数 


FX)=yE+(1 -ye)expF), y= , (10) 


= 守 Ep € {EL,6e}. (17) 
上 述 公 式 中 参数 上 的 横 杠 表示 其 为 无 量 纲 参数 . 
将 时 间 影响 因素 结合 固化 速率 és 以 及 打印 速 


度 i, 可 得 : 


X=) = 竺 -te 生地 全 (18) 


将 (18) 式 代入 (9) 式 可 得 无 量 纲 屈 曲 方程 : 


_ yal 120-vw) hn 1 
(8:W Xx) .xX 十 .ee X 二 此 “二 SxW 
hn 'ép h ” Er 
i (19) 
XWw)x=0 
éE 


w= w/h 
ws r= (a) R= (A)th, 
by = Bar = (ol) 
tt 嫩 [区 ,经 } 

(20) 


其 中 wa 为 无 量 纲 增幅 系数 ,其 目的 在 于 在 参数 固定 
时 ,利用 wu 使 得 及 = 1, 式 (19) 可 变换 为 . 

= 
dt 


(8:W XX) XK + Y=- 8 
E E E 
(21) 
其 中 大 为 无 量 纲 曲 率 ,屈曲 方程 将 被 曲率 大 与 固化 速 
率 & 控制 
k= (22) 
将 (18) 式 代入 边界 条 件 (10), 可 得 柱 壳 边界 条 件 
沿革 = -K 和 和 = 0 分别 为 : 
= es 0, i T= 0, 
i a (23) 


—c 一 —c 
wxxlx=0 =0, (8W x)xlx=0 = 0， 


ee 
w— ss(Xwx)x=0, 


2.2 ”特征 值 问题 的 求解 过 程 


效 分 方程 (21) 可 利用 伽 辽 人 金 (Galerkin) 方 法 进行 
求解 : ， 
| (Rom,)dX = 0， (24) 
其 中 余 函 数 尺 为 


一 2 
k= Kk 
vsX+12(1 -vy) gw 
E E 


R= RO) = (8,W xx) xx + 


] 
一 (XWx)X = 0. 
五 
(25) 


其 中 6%w， 为 试 函数 . Ww 可 以 设 为 满足 边界 条 件 
的 基 函 数 与 未 知 广义 坐标 C， 的 线性 组 合 , 即 wW，= 
3, CX” .在 式 (24) 中 ,每 个 广义 坐标 通过 5m，= 
(7/6C,)6C 其 得 出 的 相应 方程 必须 同时 求解 .通过 
采用 增 量 迭代 (牛顿 迭代 法 ), 对 得 到 的 广义 坐标 代数 
方程 组 进行 数值 求解 .其 为 描述 柱 壳 轴 对 称 屈曲 的 特 
征 值 问题 . 

考虑 式 (20)(22) 中 的 无 量 纲 参 数 , 柱 壳 径 厚 比 r/ 户 
厚度 及 屈曲 长 度 1,, 国 化 速率 总 , 增幅 系数 w, 以 及 
垂直 方向 打印 速度 ; 的 影响 .可 得 到 / 义 届 曲 长 度 函 
数 , /v = 1 (Ep, 


2.3” 柱 壳 的 塑性 破坏 


除 弹 性 届 曲 外 , 柱 壳 也 会 因 自 重 因素 影响 ,达到 
材料 塑性 屈服 强度 er。 而 发 生 破坏 ,其 中 的 下 标 p 代 
表 " 塑 性 破坏 "通常 来 说 ,塑性 破坏 的 关键 位 置 为 柱 壳 
底部 , 即 x = 0 处 ,此 处 因 自 重 产生 的 轴 向 应 力 达 到 最 
大 值 . 柱 壳 底部 塑性 破坏 的 屈服 极限 强度 通常 可 表示 
为 : 


pgl = |op|， (20) 


其 中 | 代表 届 服 强度 oj 的 绝对 值 ，! 为 圆柱 壳 的 长 
度 .与 弯曲 刚度 D, 变化 进行 相似 处 理 .由 于 固化 进程 
因素 影响 ,屈服 强度 在 圆柱 过 长度 增 长 方向 是 不 均匀 


的 .在 x 轴 方向 ,可 将 e 由 oo 代 蔡 ,圆柱 壳 底部 屈服 
强度 随时 间 变 化 可 以 给 定 为: 
Op(x = 0,0 =h(Dope. (27) 


其 中 oo 为 打印 进程 中 = 0 时 刻 的 屈服 强度 , = 
h(n) 为 届 服 强度 随时 间 演 化 的 固化 特征 函数 . 


由 于 塑性 届 服 极限 由 材料 所 决定 ,几何 参数 曲率 
对 其 无 影响 ,可 得 出 与 文献 2 完全 相同 的 两 种 无 量 纲 
塑性 破坏 长 度 随 固 化 速率 变化 函数 . 
(a) 线 性 固化 届 服 函数 : 
1 3 
by = 一 一 0<é,<]1, (28) 
l=~éts 
(b) 指 数 衰减 固化 届 服 函数 
i,— [ys + (1 — yo)exp(-&1,)] =0, (29) 


文献 2 利用 Newton - Raphson 迭 代 法 将 (29) 求 解 ,得 
出 一 个 封闭 的 精确 逼近 的 数值 解 ,并 将 其 拟 合 , 拟 合 


优 度 R2 介 于 0.993 与 0.999 之 间 : 
_ 元 二 了 
ly 至 p,0 1 二 和 
] 十 (EE) 
Sof (30) 
se 1.181 
Ser = rf (V5) = TO 844 
= p(ys) = 1.4660)032， 
其 中 无 量 纲 参 数 , 塑性 破坏 长 度 1,, 固化 速率 & 分别 
为 
b= 如 
(on 
gr Be 
< Or ,0 区 
c= éo el{é,é). 
pgl 
线性 固化 


届 曲 长 度 


图 6: 线性 而 化 速率 影响 下 , 柱 壳 塑性 届 服 长 度 1, 随 
固化 速率 z= 变化 曲线 图 


指数 衰减 固化 y=2 


届 曲 长 度 2 


0.0 0.5 1.0 1.5 2.0 


固化 速率 


图 7: 指数 衰减 固化 速率 影响 下 , 柱 壳 塑 性 届 服 长 度 
随 固化 速率 z= 志 变化 曲线 图 . 


如 式 (3D) 所 示 , 从 材料 性 能 和 打印 工艺 数据 中 ,可 
以 先 验 性 地 计算 无 量 纲 固化 速率 已. 如 式 (28) 所 
述 , 当 召 -1 时 ,无 量 纲 届 服 长 度 7 接近 无 穷 大 . 因 
此 ,对 于 线性 固化 过 程 ,固化 速率 趋 近 1 的 这 个 值 应 该 
被 认为 是 塑性 破坏 的 上 限 . 当 召 > 1 时 ,如 式 (3D 所 
示 ,固化 速率 &L 引起 的 届 服 强度 的 增长 快 于 打 刷 速 
度 i 所 控制 的 应 力 的 增长 .因此 ,圆柱 壳 底 部 的 应 力 不 
再 达到 届 服 强度 , 则 不 会 发 生 塑性 倒塌 . 

如 式 (30) 所 述 , 当 吾 一 co 时 其 塑性 破坏 长 度 郊 
收敛 于 二 0.y. 因此 ,对 于 指数 衰减 固化 过 程 , ,0.y 应 
该 被 认为 是 塑性 破坏 的 上 限 . 


2.4 Suiker 数 与 弹性 屈曲 与 塑性 破坏 的 Suiker 判 据 


当 塑 性 届 服 长 度 小 于 弹性 届 曲 长 度 , 即 1，< 1 
时 ,圆柱 过 会 发 生 塑性 坦 塌 ;同样 , 当 塑性 届 服 长 度 大 
于 弹性 届 曲 长 度 ， 即 凡 > 1, 时 ,圆柱 壳 会 发 生 弹 
性 届 曲 .这 一 用 于 分 析 可 能 失效 机 制 的 标准 可 以 通 
过 调用 几何 、 材 料 和 打印 参数 来 方便 地 表示 ,利用 


式 (20) 与 (31) 可 得 
os 本 
<A: elastic buckling, 
7 (32) 
>A: plastic collapse, 
ly 
二 中 
h 1 lop,ol 
= (0) es, 
0 (pg) (33) 


Le 至 1 (CE 及 ， 必 二 


Suiker2 首 次 提出 式 (33) 中 所 述 无 量 纲 参数 入 ， 
由 于 其 重要 性 , 称 A 为 Suiker 数 . 人 在 被 打印 物体 


chinaXiv:202012.00019v1 


几何 参数 ,材料 参数 固定 时 为 常数 ,不 随时 间 变 化 . 当 
评判 具体 直 壁 结构 弹性 屈曲 和 塑性 破坏 之 间 的 竞 
争 关系 时 ,完全 可 以 预先 计算 出 A, 以 便 对 破坏 机 制 
进行 预测 .在 实际 打印 过 程 中 , 常 采 用 试 错 的 方式 ,来 
对 直 壁 结构 的 失效 进行 分 析 , 即 使 在 相同 的 材料 参 
数 下 ,由 于 几何 参数 的 差异 ,失效 机 制 同样 会 发 生变 
化 .而 Suiker 数 人 的 给 出 ,为 通过 提升 材料 性 能 或 改变 
几何 参数 来 达到 材料 的 最 大 利用 率 明确 了 方向 ,这 对 
于 分 析 3D 打 印 进程 中 直 壁 结构 的 力学 性 能 及 其 失效 
机 制 方面 具有 十 分 重大 的 意义 . 


3 ” 柱 壳 失效 问题 数值 解 


在 结果 分 析 中 ,首先 通过 计算 临界 届 曲 长 度 
来 考虑 届 曲 的 情况 ,临界 届 曲 长 度 i,(Ep, 有 D 可 从 方 
程 (19) 的 弱 解 形式 与 边界 条 件 (23) 解 得 .将 代入 方 
程 弱 解 形 式 (24), 转 换 成 特征 值 问题 .所 求 最 小 特征 值 
代表 临界 届 曲 长 度 1, 通过 假设 位 移 为 多 项 式 形式 
来 进行 计算 . 


N 
h(R) = 》， CX n=1,2,.,N. (34) 


n=1 


其 中 C 代表 未 知 广义 坐标 .采用 式 (33) 给 出 的 评判 
机 制 ,可 以 进行 临界 届 曲 长 度 1,(&6, 避 与 塑性 届 服 长 
度 1,(E,) 分 别 在 线性 固化 进程 以 及 指数 衰减 固化 进 


程 中 的 竞争 关系 对 比 评判 . 


3.1 国 化 进程 影响 下 数值 求解 分 析 


对 于 圆柱 壳 在 固化 进程 影响 下 的 数值 求解 分 
析 , 通 过 调整 式 20) 增 幅 系数 w, 使 得 万 = 1 时 ,无 量 纲 
临界 届 曲 长 度 取决 于 两 个 无 量 纲 参数 ,固化 速率 如， 
j 率 天 即 筷 (局 . 其 中 泊 松 比 v = 0.3, 并 设 定 设 定 
无 量 纲 曲率 取 值 范围 为 如 <KE< +. 


线性 固化 


届 曲 长 度 7 


T T T T 1 
0.00 0.02 0.04 0.06 0.08 0.10 


固化 速率 E 


图 8， 柱 壳 在 3D 打 印 进程 中 ,由 于 线性 固化 参数 影 
响 , 届 曲 长 度 1, 随 固化 速率 E = 志 变化 曲线 图 ,其 中 
曲率 让 <k< 于 . 


由 图 8 可 以 看 出 ,在 线性 固化 进程 影响 下 , 当 曲 
率 庆 固定 时 , 届 曲 长 度 随 固化 速率 的 增加 旦 指数 形 
式 增加 ; 当 固化 速率 固定 时 , 届 曲 长 度 7, 随 曲率 的 
增 大 而 增 大 : 当 曲率 下 取 值 越 大 时 , 届 曲 长 度 随 固化 
速率 翅 的 增长 幅度 明显 增 大 ; 当 固化 速率 (& 一 0 
时 ,描述 的 为 底部 固 支 ,上 部 自由 的 均 质 柱 壳 在 自重 
作用 下 的 届 曲 行为 . 


二 


指数 衰减 固化 y=2 


屈曲 长 度 7 


T T T 1 
0.00 0.04 0.08 0.12 0.16 0. 20 


固化 速率 


图 9: 柱 壳 在 3D 打 印 进程 中 ,由 于 指数 衰减 固化 参数 
影响 , 届 曲 长 度 , 随 固化 速率 2 = 已 变化 曲线 图 , 计 
算 方便 起 见 ,其 中 曲率 专 <X< 4 y= 全 =2. 


由 图 9 可 以 看 出 ,在 指数 衰减 固化 进程 影响 下 , 当 
曲率 下 固定 时 , 届 曲 长 度 随 固化 速率 z 的 增加 而 增 
加 ,并 当 E oo 时 ,屈曲 长 度 收 敛 于 0 y; 当 固化 速 
率 (&) 固定 时 ,屈曲 长 度 i, 随 曲率 天 的 增 大 而 增 大 ; 当 
曲率 取 值 越 大 时 , 届 曲 长 度 随 固化 速率 如 的 增长 
幅度 明显 增 大 ; 当 固 化 速率 (2。 -0) 时 ,所 得 届 曲 长 
度 与 线性 固化 进程 中 完全 相同 ,同样 描述 的 为 柱 过 自 


8 


重 作用 下 的 届 曲 行为 

通过 图 7 图 8 的 分 析 可 以 看 出 ,二 者 均 描述 了 打印 
进程 中 均 质 柱 壳 到 非 均 质 柱 壳 的 自重 作用 下 届 曲 长 
度 的 变化 特性 ,以 及 曲率 对 结构 稳定 性 的 影响 , 且 符 
合 在 时 间 参 数 影 响 下 届 曲 长 度 应 随 材料 固化 速率 增 
加 而 增加 的 一 般 事实 规律 


3.2 ”弹性 屈曲 与 塑性 破坏 数值 求解 比较 分 析 


对 于 柱 壳 弹性 屈曲 与 塑性 破坏 数值 求解 ,分 析 线 
性 固化 进程 与 指数 衰减 固化 进程 下 的 弹性 届 曲 与 塑 
性 破坏 竞争 关系 。 

线性 固化 进程 情况 下 ,选取 曲率 k = 证 ,并 对 此 
参数 所 求 数值 解 进行 拟 合 ,可 用 近似 函数 表示 ,其 拟 
合 优 度 R2 = 0.999 


i,, = 8.7192 + 3.0652 exp(18.2639B)， (35) 


其 中 当 恕 趋 近 于 0 时 , 即 无 量 纲 临界 屈曲 长 度 与 国 化 


吧 | 


塑性 地 塌 ,反之 将 发 生 弹 性 届 曲 .弹性 屈曲 与 塑性 南 
塌 的 竞争 关系 明显 依赖 于 B 值 的 变化 . 当 p 取 值 越 大 
时 , 标 度 比 7/1, 曲线 下 降 越 快 .对 于 图 9 而 言 , 当 B > 6 
时 ,只 要 入 > Lo = 11.65 时 , 塑性 拥 塌 将 不 会 发 生 , 弹 
性 届 曲 为 唯一 的 失效 机 制 |. 
当 A < lo = 11.65 时 ,塑性 破坏 将 发 生 在 当 
固化 速率 较 小 时 ,弹性 届 曲 将 发 生 在 固化 速率 二 
较 大 时 . 随 着 束 的 增 大 ,1,/1, 向 x 负 轴 方向 减 小 , 直 
到 万/ 万 = 0， 弹性 届 曲 为 唯一 失效 机 制 .具体 来 
说 , 当 B = 5 时 , 当 A > (yl)wax = 14.58 时 ， 弹 性 
屈曲 为 唯一 失效 机 制 。 当 0 < A < (fly)mar 时 ，0 < 
如 < 4b < 如 为 弹性 届 曲 , a < 你 < 上 5 为 逆 性 破坏 , a 
b 可 由 p(y/iy) 求 得 , 其 中 (4,/i,) = 人 

指数 衰减 固化 进程 情况 下 ,对 6, i 所 求 数值 解 
进行 拟 合 ,可 用 近似 函数 表示 ,其 拟 合 优 度 R?2 = 0.999. 


一 


1 = 22.9168 - 11.4916 exp(—8.0559e7). (37) 


其 中 当 况 = 0 即 无 量 纲 临界 届 曲 长 度 与 固化 速率 


拟 合 近 似 函 数 ,对 弹性 届 曲 与 塑性 坊 塌 的 竞争 关系 进 
行 表述 . 


Re] 
"|S 


* = (8.7192 + 3.0652 exp(18.26392))(1 ~ Be), (36) 
p 


其 中 p = 如 有 反 . 


站 


| 线性 固化 


Som 


标 度 比 和 /7 


T T 
0. 00 0.05 0.10 0.15 0.20 


图 10: 线性 固化 进程 中 ,圆柱 过 打印 失效 图 .长 度 标 
度 比 /i, 在 选 定 的 固化 比 下 , 随 固化 速率 2 = 十 变 
化 曲线 图 .其 中 p = 如 / 束 , 其 范围 从 5 到 10. 此 曲线 由 


无 关 时 ,16.o = 11.65. 同时 可 用 拟 合 近似 函数 ,结合 
式 (30) (37) 对 弹性 屈曲 与 塑性 破坏 的 竞争 关系 进行 表 
述 ,可 得 


Jer 一 e 
= (22.9168 — 11.4916 exp(—8.055985)) 
p 
(38) 
S| 
1 和 = 
1+ EE) 
Eef 
25 一 一 站 .1 
pr 0.2 
指数 衰减 固化 y=ys=y。=2 Es 0 
20 E 
Sa st 10 
-== 20 
| -一 40 
15 


o 1 1 
固化 速率 


图 11: 指数 衰减 固化 进程 中 ,圆柱 过 打印 失效 图 .长 度 


式 (3 人 绘制 ,其 中 入 代表 了 线性 固化 进程 中 ,弹性 屈曲 
与 塑性 破坏 之 间 的 竞争 关系 . 


如 图 10 所 示 , 其 中 6 的 取 值 为 ,5,6,7,8,9,10. 具 体 来 
说 , 当 标 度 比 iy1, 曲线 上 的 点 大 于 人 的 值 时 ,将 发 生 


标 度 比 1,/i, 在 选 定 的 固化 比 下 , 随 固化 速率 2 = 普 
变化 曲线 图 .其 中 p = 名 /名 ,其 范围 从 0.1 到 40. 此 曲线 
由 式 (38) 绘制 ,其 中 大 代表 了 指数 衰减 固化 进程 中 , 弹 
性 届 曲 与 塑性 破坏 之 间 的 竞争 关系 . 
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如 图 10 所 示 , x 坐标 采用 对 数 坐 标 .其 中 p 的 取 
值 为 0.1，0.2，0.5，1，2，5，10，20，40. 具 体 分 析 
与 上 述 类 似 , 当 标 度 比 L,y1, 曲线 上 的 点 大 于 入 的 值 
时 .将 发 生 塑性 破坏 ,反之 将 发 生 弹 性 屈曲 .如 图 10 所 
示 , 当 6 > 20 时 ,只 要 A > 11.65, 塑性 破坏 将 不 会 
发 生 . 弹 性 屈曲 为 唯一 的 失效 机 制 , 当 A < 11.65 时 ， 
0 < 可 < wp < 如 为 塑性 破坏 ,a < 或 < 5 为 弹性 届 
曲 , a, 5 可 由 各 (efiy) 求 得 


具体 来 说 , 当 B = 0.1 , 当 A > (6/1y)mar = 22.54 
时 ,弹性 届 曲 为 唯一 失效 机 制 . 当 0 < A < (oy/ly)max 
时 ,0 < Eh < a,b < 4 为 弹性 届 曲 , a < E6 <5 为 塑性 
破坏 , a, b 可 由 训 (1/iy) 求 得 ,其 中 4/1,) = 
通过 图 9 图 10 的 分 析 可 以 看 出 ,固化 比 8 对 弹性 
届 曲 与 塑性 破坏 之 间 的 竞争 关系 影响 巨大 . 当 p 越 小 


4.1 弹性 屈曲 


图 12: 柱 壳 届 曲 生长 示意 图 


图 11 中 描绘 了 柱 这 由 于 在 X 轴 方 向 以 恒定 速 
率 生 长 ,结构 从 最 初 小 变形 ， 其 次 达到 临界 届 曲 高 
度 而 发 生 失 稳 直 至 结构 失效 的 示意 图 .针对 柱 壳 结 
构 弹 性 屈曲 分 析 中 所 使 用 的 参数 列 于 表 1, 对 在 自 
重 作用 下 屈曲 的 力学 行为 进行 析 , 其 中 选择 了 线性 
固化 以 及 指数 衰减 固化 两 种 常见 的 固化 方式 进行 


时 , 越 容易 发 生 塑 性 破坏 ; 当 pB 越 大 时 , 越 容易 发 生 弹 
性 屈曲 .而 标 度 率 曲 线 与 A 值 相交 处 ， 代 表 两 种 破 
坏 机 制 的 临界 位 置 .通过 对 标 度 比 曲 线 的 变化 趋势 发 
现 ,线性 固化 进程 中 , 当 B = 5 时 , 标 度 比 1/1; 为 先 增 
加 在 减少 ;而 当 6 > 5 时 , 标 度 比 j/ 为 单调 递减 .在 
间 数 衰减 固化 进程 中 , 当 B < 5 时 , 标 度 比 心 /7 为 先 
增加 再 减少 ; 当 B > 10 时 , 标 度 比 i,/i; 为 先 减少 再 增 
加 .由 此 可 以 看 出 ,固化 比 6 对 标 度 比 曲线 的 变化 趋势 
也 有 很 大 的 影响 , 即 在 给 定 几何 参数 下 ,固化 比 6 决定 
了 弹性 屈曲 与 塑性 破坏 这 两 种 失效 机 制 的 主导 地 位 . 


人 


4 给 定 打印 参数 下 的 影响 分 析 


本 节 首 先 使 用 商用 有 限 元 软件 ABAQUS 对 3D 打 
印 进程 中 柱 壳 结构 的 届 曲 行为 进行 了 模拟 ,对 上 
文 所 得 屈曲 参数 模型 进行 了 验证 .随后 采用 了 
与 Suikert24 完 全 相同 的 打印 参数 ,将 圆柱 党 结构 与 其 
所 研究 的 b = 625mm,d = 250mm 矩形 布置 直 壁 结构 
进行 对 比分 析 . 最 后 针对 该 类 混凝土 所 打印 圆柱 壳 结 
构 的 失效 机 制 , 即 弹 性 屈曲 与 塑性 破坏 之 间 的 竞争 关 
系 进行 了 有 具体 分 析 . 


模拟 .根据 式 (14)(20)， 对 线性 固化 进程 选择 天 全 = 
0.02，0.05，0.1, 对 指数 衰减 固化 进程 选择 了 总 = 
0.02, 0.1, 0.2 进行 模拟 ,以 便 对 参数 模型 进行 验证 . 


表 1 Suiker 打 印 固化 参数 


参数 值 
壳 体 厚度 h = 43.5(mm) 
单 层 高 度 11 = 9.2(mm) 
材料 密度 p = 2020(kg/m3) 
喷头 速度 yn = 83.3(mm/s) 
泊 松 比 y = 0.3(-) 
线性 固化 
初始 弹性 模 量 Eo = 78100(Pa) 
固化 速率 天 = 0.02,0.05,0.1(-) 
指数 衰减 固化 
初始 弹性 模 量 Eo = 78100(Pa) 
固化 速率 下 = 0.02,0.1,0.2(-) 


注意 ,为 贴 合 实际 打印 尺寸 以 及 数值 计算 方 
便 , 式 (20) 中 的 c 为 无 量 纲 系数 , 当 实际 打印 壳 体 厚度 
为 43.5mm 时 ,通过 调节 au = 8.297 来 确保 无 量 纲 壳 体 
厚度 h=1. 其 中 ,在 参数 一 定 的 情况 下 ,无 量 纲 固化 速 
率 E. 与 固化 速率 &j 为 正 线 性 关系 . 即 无 量 纲 固化 速 
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率 &j; 愈 大 ,固化 速率 és 也 愈 大 .由 式 (33)4 可 得 , 当 无 
纲 壳 体 厚 度 h = 1 时 ,无 量 纲 屈曲 长 度 心 由 无 量 纲 
化 速率 &s, 以 及 无 量 纲 曲率 大 决定 , 即 16, = cr 人 er 月. 
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13: 线性 固化 参数 影响 下 ,FEM 结 果 与 数值 计算 结 
果 对 比 图 


由 图 12 可 以 看 出 :在 线性 固化 参数 影响 下 ,得 出 
的 有 限 元 分 又 屈曲 模拟 与 参数 模型 的 结果 整体 较为 
接近 .但 是 随 着 总 的 增 大 ,有 限 元 结果 与 参数 模型 结 
果 差距 增 大 .具体 来 说 , 当 吏 = 0.02 时 ,二 者 相差 3.1%; 
当 翅 = 0.05 时 ,二 者 相差 8.6%; 当 翅 = 0.1 时 ,二 者 相 
差 25%. 
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14: 指数 衰减 固化 参数 影响 下 ,FEM 结 果 与 数值 计 
算 结 果 对 比 图 


由 图 13 可 以 看 出 :在 指数 衰减 固化 参数 影响 
下 ,对 于 有 限 元 模拟 ,得 出 的 有 限 元 模拟 分 又 届 曲 结 
果 与 参数 模型 的 结果 整体 同样 较为 接近 .但 是 随 着 
的 增 大 ,有 限 元 结果 与 参数 模型 结果 差距 增 大 .具体 
来 说 , 当 至 = 0.02 时 ,二 者 相差 1.9%; 当天 = 0.1 时 ,二 
者 相差 7.2%; 当 瑟 = 0.2 时 ,二 者 相差 10%. 

当 = 625mm,d = 250mm 时 ,曲率 应 为 
224D ~ 1/6.4， 其 中 柱 壳 厚度 为 43.5mm， 单 
打印 高 度 为 9.2mm， 根 据 式 (11)(14)(40) 可 得 如 
0.000256(s-1),i = 0.000438(m/s)， 代入 式 (20)4 可 得 ， 
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Ej = 0.0254. 而 对 于 此 类 混凝土 , 当 曲 率 4 < 
有 < 证 时 ,无 量 纲 固化 速率 你 的 范围 由 0.0199 增长 
至 0.0397, 在 此 范围 内 ,有 限 元 模拟 分 又 届 曲 结果 与 参 
数 模型 结果 相差 不 大 ,说 明 参 数 模型 能 对 混凝土 3D 打 
印 柱 壳 的 分 又 届 曲 行为 进行 较 好 的 模拟 

在 经 典 理论 中 ,对 于 柱 壳 屈曲 的 控制 方程 的 解 
采用 正弦 近似 的 方式 写 出 级 数 形式 : w = 2 YA4. 
sin(ax).sin(By), 其 解 沿 圆 周 方向 同样 采用 屈曲 波形 为 
正弦 形式 ,此 解 更 加 符合 实际 情况 .而 本 文 在 第 一 节 
的 理论 推导 中 ,考虑 到 混凝土 3D 打印 圆柱 壳 由 于 殿 
初始 抗 弯 刚度 Eo 较 低 ,不 可 能 在 较 大 挠 度 下 仍 保持 
完好 ,所 以 构建 的 为 轴 对 称 情况 下 的 参数 模型 ,将 解 
设置 为 w = w(x). 在 对 于 总 值 较 大 时 ,混凝土 3D 打 印 
柱 壳 届 曲 长 度 更 高 ,而 距离 壳 体 底部 距离 越 大 处 ,其 
水 平 届 曲 形态 的 差异 越 大 ,有 限 元 模拟 分 叉 届 曲 结 
与 参数 模型 结果 相差 也 越 大 . 


4.2 对比 分析 


根据 Suiker22 研 究 发 现 , = 625mm,d = 250mm 
矩形 布置 直 壁 结构 ,其 失效 机 制 的 上 界 由 塑性 破坏 
决定 (而 不 是 侧面 固 支 下 的 弹性 屈曲 ), 下 界 由 侧面 简 
文 下 的 弹性 屈曲 决定 ,具体 参数 列 于 表 1, 固 化 进程 
式 (39)(40) 决 定 . 当 采 用 了 与 Suiker 完 全 相同 的 参数 
时 ,曲率 上 = 1/6.4, 即 有 hh = 43.5mm, r = 278.4mm. 

Suiker 通 过 对 直径 为 70mm， 高 度 为 140mm 
的 4 到 6 个 圆柱 形 混凝土 试 件 进行 单 轴 压缩 试验 , 测 
量 出 混凝土 弹性 刚度 与 抗 压 屈服 强度 在 不 同时 间 点 
的 实验 数据 ,对 此 实验 数据 采用 最 小 二 乘法 进行 拟 
合 ,得 到 
(1) 弹 性 刚度 (Pa) 随 固化 时 间 z(s) 线性 变化 函数 , 拟 
合 优 度 R2 = 0.995. 


BE,(t) = 78100 + 201, (39) 


(2) 抗 压 届 服 强度 全,,(Pa) 随 固化 时 间 xs) 线性 变化 函 
数 , 拟 合 优 度 R2 = 0.991. 


Spx(t) = 5984 + 2.451. (40) 


在 此 打印 进程 参数 与 几何 参数 影响 下 ,侧面 简 支 的 矩 
形 布 置 直 壁 结构 屈曲 高 度 = 0.210m, 侧面 固 支 情 
况 下 人 = 0.431m, 塑性 破坏 长 度 1, = 0.421m. 因此 矩 
形 布置 的 直 壁 结构 实际 失效 尺寸 范围 应 为 0.210m < 


1 < 0.421m. 进一步 由 纵横 比 yd = 2.5 确定 出 侧 壁 
对 主 壁 提供 的 转动 刚度 ,从 而 得 出 更 加 确切 的 屈曲 高 
度 16 = 0.269m. 
在 参数 完全 相同 情况 下 ,混凝土 3D 打 印 柱 壳 的 
无 量 纲 屈曲 高 度 /， = 13.3 等 同 于 1 = 0.578m， 由 
式 (31) (28) 可 以 得 出 ,此 时 如 = 0.282, 圆柱 壳 无 量 纲 
塑性 破坏 长 度 1, = 1.393 等 同 于 1, = 0.421. 即 当 曲 
率 k x 1/6.4 时 , 柱 壳 结构 会 因为 达到 抗 压 届 服 强 度 
而 发 生 塑性 破坏 . 
令 柱 壳 届 曲 高 度 与 Suiker 计 算 相 同时 , 即 1， = 
0.269m 等 同 于 1 = 6.184, 带 入 式 (43), 可 得 曲率 = 
17.74. 即 代表 着 相同 的 稳定 性 要 求 时 , 柱 壳 结构 能 比 
矩形 布置 直 墙 结构 拥有 更 大 的 框架 . 
显然 , 相 比 于 和 矩形 布置 的 直 壁 结构 ,由 于 结构 的 
完整 性 (相当 于 侧面 完全 固 支 ) 以 及 曲率 的 引入 , 柱 壳 
结构 能 更 加 充分 利用 材料 , 即 在 单 层 打印 材料 用 量 相 
同 的 同时 , 柱 壳 结构 能 达到 更 高 的 层 数 ,具有 更 强 的 
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43 失效 机 制 分 析 


具体 来 说 ,对 于 柱 壳 限定 的 几何 参数 , 当 h = 
43.5mm 时 ,为 了 计算 方便 ,通过 调整 a 的 值 ,使 得 h = 
1, 通过 式 (20) 计 算 可 得 , ac = 8.297, 即 (a 铁 )* = 22.99. 
由 此 可 得 76, = 22.99 1,. 

对 于 限定 的 混凝土 材料 及 固化 参数 , 当 其 弹性 刚 
度 和 抗 压 届 服 强 度 由 式 (39)(40) 给 定时 , 即 此 类 参数 
均 被 给 定 , 如 = 对 和 ,= 0.000256(s-),é6 = 只 = 
0.00041(s-!),p = 2020(kg/n5),g = 9.81lm/s*. 由 
式 (20)(31) 可 以 看 出 ,唯一 决定 无 量 纲 固化 速率 & 的 只 
9 柱 壳 轴 疝 增长 速率 i. 打印 材料 从 喷头 每 单位 时 间 
排出 的 体积 2 = vihu, 每 层 材 料 打 印 所 需 时 间 7 = 
2xr/vn, 结合 式 (11), 可 得 


0 tvyn hh h 
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结合 式 (41), 并 带 入 上 述 参 量 可 得 出 无 量 纲 固 化 速 
率 #4 与 &。 随 曲率 k 的 变化 函数 


-| 1 Do 1iég 


= (一 3 2 = 0.004k, 
6 则 
5 和 jcaono ) 
< = 一 =0.04409K， 
pgl 


所 得 散 点 图 ,进行 指数 拟 合 , 拟 合 优 度 R2 = 0.999 


1 = 5.63943 + 44.07928 exp(—0.285097). (43) 


利用 式 (33), 可 计算 出 此 类 参数 下 ,Suiker 数 和 人 = 
6.94196. 结合 式 (28)(42), 可 得 无 量 纲 塑性 破坏 长 度 i， 
随 无 量 纲 半径 7 变化 函数 .将 1,() .入 与 式 (43) 绘 制 于 
图 14. 
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图 15: 圆柱 党 混凝土 3D 打 印 ,在 混凝土 参数 固定 , 打 
印 参数 固定 的 影响 下 ,塑性 破坏 长 度 1, 与 弹性 屈曲 长 
度 心 随 无 量 纲 半径 7 变化 图 ， 


如 图 12 所 示 , 当 打印 参数 由 式 (39)(40) 与 表 1 决 定 
时 ,Suiker 数 和 也 随 之 固定 .而 i(7) .入 与 i,(7) 两 条 曲 
线 的 交点 4， 代 表 着 塑性 破坏 与 弹性 屈曲 两 种 破坏 
机 制 的 竞争 关系 的 临界 点 , 即 此 时 7 = 1,. 由 图 12, 可 
以 看 出 两 条 曲线 的 交点 处 , 径 厚 比 7 = 7.73496, 失效 
长 度 1 = 10.53515, 其 等 同 于 曲率 = 1/7.73496, 失 
效 长 度 0 = 0.4582m. 当 无 量 纲 半径 7 < 7.73496 时 , 塑 
性 破坏 为 结构 的 失效 机 制 , 当 7 > 7.73496 时 ,弹性 届 


将 式 (42)1, 以 步 长 为 0.5, 曲 率 k 范围 从 1/5 到 1/10, 计 多 


曲 为 结 构 的 外 效 机 制 |. 
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图 16: 圆柱 壳 混 凝 土 3D 打 印 ,在 混凝土 参数 固定 , 打 
印 参数 固定 的 影响 下 ,有 效 长 度 与 失效 长 度 区 域 区 分 
图 . 


如 图 13 所 示 , 当 弹性 屈曲 长 度 心 与 塑性 破坏 长 
度 岂 通过 Suiker 数 人 联系 起 来 时 ,可 将 有 效 长 度 与 失 
效 长 度 二 者 区 域 进 行 区 分 ,弹性 屈曲 与 塑性 破坏 二 
者 区 域 进行 区 分 .图 13 中 , 蓝 色 区 域 B 为 有 效 长 度 区 
,黄色 区 域 C 为 塑性 破坏 区 间 , 红 色 区 域 D 为 弹性 
届 曲 区 间 . 在 参数 固定 的 实际 柱 壳 打印 过 程 中 ,通过 
选择 曲率 大 小 ,或 者 通过 在 不 同 失效 区 域 采 用 相对 应 
的 辅助 手段 进行 结构 增强 的 方式 ,以 达到 预期 打印 高 
度 
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5 结 论 


为 了 提高 3D 打 印 直立 墙 的 整体 稳定 性 ,我 们 提 
出 了 一 个 简单 可 行 的 普 适 方案 ,， 即 适当 增加 直立 墙 
治 轮廓 方向 的 曲率 将 直立 墙 沿 轮 廓 方向 打印 成 波纹 
形状 。 直 立波 纹 墙 每 个 波纹 可 以 是 圆柱 形 也 可 以 不 
是 , 同时 其 曲率 也 可 以 变化 。 为 了 验证 我 们 的 方案 ， 
本 文 在 Suiker 等 人 工作 的 基础 上 ， 通 过 增加 曲率 把 
直立 板 墙 扩 展 到 直立 波纹 墙 ， 并 对 直立 柱 壳 墙 整体 
稳定 性 进行 了 细致 研究 。 发 现 柱 壳 形 直立 墙 的 临界 
届 曲 高 度 是 矩形 直 墙 的 2 倍 多 , 表明 增加 直立 墙 沿 轮 
廓 方向 的 曲率 是 提高 直立 墙 整体 稳定 性 的 合理 方案 。 
力学 上 ， 这 个 方案 之 所 以 可 以 提高 直立 墙 整体 稳定 
性 本 质 上 是 通过 增加 曲率 提高 了 直立 墙 的 整体 拓扑 
刚度 ， 是 提高 直立 墙 整体 稳定 性 的 普 适 方法 ， 而 直 
立波 纹 墙 结构 是 提高 直立 墙 整体 稳定 性 的 简单 可 行 
的 构 形 。 具体 可 总 结 出 以 下 六 个 结论 
(1) 参 数 模 型 描述 的 是 打印 进程 中 均 质 圆柱 过 到 非 均 
质 圆柱 这 在 自重 作用 下 届 曲 长 度 的 变化 特性 , 即 用 


公式 描述 了 圆柱 壳 逐 层 打印 过 程 中 ,结构 力学 性 能 
在 空间 中 不 均匀 的 情况 下 ,得 出 其 屈曲 行为 的 分 析 结 
果 , 且 符合 一 般 事实 规律 . 
(2) 在 具体 研究 中 ,本 文 将 每 个 波纹 视 为 单 胞 进行 和 
究 , 取 单 胞 为 柱 壳 形 直立 墙 时 发 现 ,其 屈曲 长 度 对 于 
率 上 有 着 较 大 的 敏感 性 , 随 着 曲率 的 减 小 ,相同 的 无 
量 纲 固化 速率 增 量 dé, 将 会 导致 更 大 的 屈曲 长 度 增 
增 量 dl. 
(3) 标 度 比 曲线 图 (9)(10) 描 述 的 为 标 度 比 i,yis 曲线 随 
固化 速率 上 的 变化 趋势 ,可 用 于 研究 分 析 两 种 破坏 机 
制 的 竞争 关系 .由 式 (33)3 所 决定 的 Suiker 数 入 描写 的 
是 塑性 破坏 与 弹性 屈曲 分 界 ,当心 /2 > A, 为 塑性 破 
坏 , 反 之 则 为 弹性 屈曲 .其 中 国 化 比 6 参数 对 标 度 比 曲 
线 的 变化 趋势 有 很 大 的 影响 .其 中 和 A 由 材料 的 初始 参 
数 决 定 ,固化 比 6 由 材料 的 固化 参数 决定 .这 两 个 参数 
决定 了 弹性 屈曲 与 塑性 破坏 这 两 种 失效 机 制 的 主导 
地 位 . 
(4) 在 数值 结果 与 有 限 元 结果 的 对 比分 析 结 果 的 差 
异 ,我 把 其 归结 于 我 所 建立 的 参数 模型 的 简化 ,对 于 
真实 混凝土 打印 参数 ,还 是 能 较 好 的 符合 . 
(5) 采 取 了 定量 分 析 的 方式 ,加 入 曲率 大 的 影响 ,对 比 
完了 文献 2 中 珑 形 布置 直 墙 结构 与 柱 壳 结 构 在 3D 打 
印 进程 中 的 失效 形式 ,发 现 柱 壳 比 矩形 布置 直 墙 结 
构 的 稳定 性 提升 了 2.15 倍 .与 此 同时 ,如 图 (13) 所 示 , 利 
用 Suiker 数 入 , 划分 出 了 在 一 定 曲 率 范 围 内 的 有 效 长 
度 与 失效 长 度 区 域 区 分 图 ,可 针对 性 的 进行 结构 增 
强 
(6) 如 图 (14) 所 示 , 利 用 Suiker 数 入 , 可 划分 出 一 定 曲率 
范围 内 的 有 效 长 度 与 失效 长 度 区 域 区 分 图 ,在 此 基础 
上 ,可 针对 性 的 对 结构 进行 增强 . 
在 未 来 的 工作 中 ,将 圆柱 形 单 胞 结构 扩展 至 变 曲 
率 以 及 整个 波纹 墙 进 行 研究 .针对 波纹 墙 中 每 个 波 
纹 ,在 参数 给 定 的 情况 下 ,理论 上 存在 材料 使 用 率 最 
高 ,结构 稳定 性 最 强 的 一 个 波形 最 优 解 . 


参 者 文献 


1 Paul G M , Rezaienia A , Wen P , et al. Medical Applications 
for 3D Printing: Recent Developments. Missouri Medicine, 2018， 
115(1):75-81 

2 Mathias W , Sebastian T, Christoph H . Rapid manufacturing of au- 
tomotive polymer series parts: A systematic review of processes, 
materials and challenges. Additive Manufacturing, 2020, 36 

3 Wong KV ,Hernandez A .A Review of Additive Manufacturing. 
lIsrn Mechanical Engineering, 2012, 2012:30 - 38 


13 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


2 


2 


4 Kumar LJ, Krishnadas Nair C G . Current Trends of Additive Man- 
ufacturing in the Aerospace Industry. Advances in 3D Printing & 
Additive Manufacturing Technologies , 2017:39-54 

Le-Bail A , Maniglia B C , Le-Bail P . Recent advances and future 
perspective in additive manufacturing of foods based on 3D printing. 
Current Opinion in Food ence , 2020, 35:54-64 


6 Roussel N,Spangenberg J, Wallevik J ,etal. Numerical simulations 


un 


of concrete processing: From standard formative casting to additive 

manufacturing. Cement and Concrete Research, 2020, 135:106075 

BosF, Wolfs R , AhmedZ ,et al. Additive manufacturing of con- 

crete in construction: potentials and challenges of 3D concrete print- 

ing. Virtual Physical Prototyping, 2016:209-225 

Astm standard f2792, standard terminology for additive manufactur- 

ing technologies. 2013. URL 

Bourell D , Chen Y , Zhou C , et al. A layerless additive manu- 

facturing process based on CNC accumulation. Rapid Prototyping 

Journal, 2011, 17(3):218-227(10) 

0 Williams C B , Mistree F , Rosen D W. A Functional Classification 
Framework for the Conceptual Design of Layered Manufacturing 
TechnologiesW ASME 2008 International Design Engineering Tech- 


~1 


Oo 


Ro 


nical Conferences and Computers and Information in Engineering 
Conference. 2008 

1 Gao W ,Zhang YB , Devarajan R,The status, challenges, and future 
of additive manufacturing in engineering. Computer-Aided Design, 
2015, 69:65-89 

2 David L B, Ming CL ,and David W R . Roadmap for Additive 
Manufacturing, University of Texas at Austin, Austin TX, 2009 

3 Fatemeh H, Farhad A . Additive manufacturing of cementitious 
composites: Materials, methods, potentials, and challenges. Con- 
struction and Building Materials, 2019, 218:582-609 

4 陈 玲 钰 ， 张 朝 弹 . 基于 BIM 技 术 的 3D 打 印 装配 式 建 筑 应 

7 


上 探讨 //2020 国 际 绿色 建筑 与 建筑 节能 大 会 论文 集 ,第 十 

国际 绿色 建筑 与 建筑 节能 大 会 暨 新 技术 与 产品 博览 会 ,苏州 

市 ,2020 年 8 月 .北京 : 中 国 城市 出 版 社 ,2020. 659-661 

5 王 香港 , 王 申 , 贾 鲁 涛 等 . 3D 打 印 混 凝 土 技 术 在 新 冠 肺炎 防疫 方 

舱 中 的 应 用 . 混凝土 与 水 泥 制 品 , 2020, 000(004):1-4,13 

6 陈 董 BIM+3D 打 印 的 装配 式 建筑 技术 组 合 应 用 的 优 
势 . 武汉 交通 职业 学 院 学 报 , 2020, 022(001):81-84 

7 侯 泽 宇 , 张 宇 , 张 超 等 . 3D 打 印 混凝土 力学 性 能 试验 方法 . 混凝土 
与 水 泥 制 品 , 2019, 283(11):5-9 

8 雷 波 , 解 文 峰 . 一 种 3D 打 印 混凝土 材料 的 试验 研究 . 新 型 建筑 材 
料 , 2020, 047(002):72-76 


9 马 国 伟 , 柴 艳 龙 , 王 里 等 .3D 打 印 陶 砂 轻 质 混凝土 的 制备 与 力学 性 
能 测试 .实验 力学 ,2020,35(1):58-66 

0 马 国 伟 , 柴 艳 龙 , 王 里 等 .3D 打 印 陶 砂 轻 质 混凝土 的 制备 与 力学 性 
能 测试 .实验 力学 ,2020,35(1):58-66 

1 段 严 , 秦 先 涛 .3D 打 印 混凝土 相关 性 能 研究 进展 .混凝土 与 水 泥 制 


品 ,2020(09):5-10 


22 


23 


24 


25 


26 


27 


28 


29 


30 


31 


32 


33 


34 


33 


36 


37 


38 


39 


A.S.J. Suiker. Mechanical performance of wall structures in 3D 
printing processes: Theory, design tools and experiments. Interna- 
fional Journal of Mechanical Sciences, 2018, 137:145-170 

Bos F, Wolfs R R., Ahmed Z Z, et al. Additive manufacturing of 
concrete in construction: potentials and challenges of 3D concrete 
printing. Virtual and Physical Prototyping, 2016, 11(3): 209-225 

V. Mechtcherine, F.P. Bos, A. Perrot, et al.Extrusion-based addi- 
tive manufacturing with cement-based materials - Production steps, 
processes, and their underlying physics: A review. Cement and Con- 
crete Research, 2020, 132 

Timoshenko SP. Strength of Materials. Lancaster Press, Inc., Lan- 
caster, PA; 1947 

Wolfs RJ M,BosFP, Salet TA M .Early age mechanical be- 
haviour of 3D printed concrete: Numerical modelling and experi- 
mental testing. Cement & Concrete Research , 2018, 106:103-116 
Vardoulakis I, Sulem J. Bifurcation Analysis in Geomechanics. 
Black Academic Professional, London; 1995 

Lim C w, Ma Y F. Computational p-element method on the effects 
of thickness and length on self-weight buckling of thin cylindrical 
shells via various shell theories. Computational Mechanics , 2003, 
31:400 - 408 

Suiker A SJ, Wolfs RJM ,Lucas S M ,etal. Elastic buckling and 
Plastic collapse during 3D concrete printing. Cement and Concrete 
Research, 2020, 135:106016 

Wolfs RJM, Suiker A SJ. Structural failure during extrusion-based 
3D printing processes. International Journal of Advanced Manufac- 
turing Technology, 2019, 104(1-4):1-20 

视 恩 淳 ,Mandal P , Calladine CR . 轴 压 圆柱 薄 壳 的 届 
木工 程 学 报 , 2001(03):20-24+30 

ReddyJN , 工 JHS . General buckling of stiffened circular cylindri- 
cal shells according to a layerwise theory. Computers & Structures, 
1993, 49(4):605-616 

周 承 个 . 弹性 稳定 理论 . 四 川 和 人民 出 版 社 , 1982 

Timshenko S P, Gere J M. Theory of Elastic Stability.2nd ed. New 
York: Dover Publications Inc,2009 

Timoshenko SP Woinowsky-Krieger S$. Theory of Plates and Shells. 
McGraw-Hill Book Company, Singapore, 1959 

Calladine C R , Barber J N . Simple Experiments on Self-Weight 
Buckling of Open Cylindrical Shells. Journal of Applied Mechanic- 
s, 1970, 37(4):1150 

Miller J R, Weaver P M, JOHNS DJ. Self-weight-buckling of verti- 
cal circular cylindrical shells. AIAA Journal, 1973 

Kumar, Yajuvindar Lal,R. The vibration and buckling of freely sup- 


分 析 . 土 


比 


ported non-homogeneous orthotropic conical shells subjected to d- 
ifferent uniform pressures. Journal of Applied Mechanics,201]1, 
061012 

Wang Y Q ,Liu YE. Free vibration and buckling of polymeric 
shells reinforced with 3D graphene foams. Results in Physics, 2019, 


14:102510 


